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Abstract. The present paper is related to a talk presented during the Symposium 
on Coherent Control and Ultracold Chemistry held during the Sixth Congress of 
the International Society for Theoretical Chemical Physics (ISTCP-VI, July 2008). 
The talk was entitled "Electronic structure properties of alkali dimers and trimers. 
Prospects for alignment of ultracold molecules". Here we report on the electrostatic 
repulsion forces of the ionic cores at short separation, involved when the potential 
energy surfaces of alkali trimers are calculated with a quantum chemistry approach 
based on effective large-core potentials for ionic core description. We demonstrate 
that such forces in the triatomic molecule can be obtained as the sum of three pairwise 
terms. We illustrate our results on the lowest electronic states of Csa, which are 
computed for the first time within a full configuration interaction based on a large 
Gaussian basis set. As a preliminary section, we also propose a brief introduction 
about the importance of alkali trimer systems in the context of cold and ultracold 
molecules. 
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1. Introduction: a brief overview on cold molecules 

The research field on translationally cold (~ IK and lower) and ultracold (~ ImK and 
lower) molecules is continuously expanding in many directions, involving an increasing 
number of groups throughout the world. The availability of gaseous samples of cold 
neutral molecules opens entirely new avenues for fascinating researches, offering the 
possibility to control all degrees of freedom of a quantum system. Molecular ions can 
also be trapped and cooled down inside atomic ion traps behaving like ionic crystals, 
where each ion is kept at a specific position with a residual motion equivalent to a 
temperature smaller than 1 K [1]. Cold neutral molecules brought new perspectives 
in high- resolution molecular spectroscopy [21 [3l HI [5]. The expected accuracy of the 
envisioned measurements with ultracold molecules makes them appear as a promising 
class of quantum systems for precision measurements related to fundamental issues: 
the existence of the permanent electric dipole moment of the electron [6], 13 El IS] 
related to CP-parity violation [TOl HU [T2] , and the time-independence of the electron- 
to- nuclear and nuclear-to- nuclear mass ratios [131 [IH [13 IISI [13 [H], or of the fine- 
structure constant [19]. Various proposals have been suggested for achieving quantum 
information devices [201 [211 [221 123] based on cold polar molecules (i.e. molecules 
exhibiting a permanent electric dipole moment). Elementary chemical reactions at very 
low temperatures could be manipulated by external electric or magnetic fields, therefore 
offering an extra flexibility for their control [23] . The large anisotropic interaction 
between cold polar molecules is also expected to give rise to quantum magnetism 
[25], and to novel quantum phases [261 [22] ■ The achievement of quantum degeneracy 
with cold molecular gases [281 1291 [301 [Ml [321 [331 [32 [35] together with the mastering 
of optical lattices [36] built a fantastic bridge between condensed matter and dilute 
matter physics. Indeed, the smooth crossover between the Bose-Einstein condensation 
(BEC) of fermionic atomic pairs and the Bardeen-Cooper-Schrieffer (BCS) delocalized 
pairing of fermions related to superconductivity and superfluidity, has been observed 
experimentally [35l [371|38l [39] . 

Cold and ultracold neutral molecules can be formed along two paths: 

• Existing molecules in their vibrational ground state can be slowed down and cooled 
by various means, such as interactions with external electric and magnetic flelds, 
or by collision with surrounding particles (see for instance the review articles in 
refs. [lOl EH Hi]. This approach concerns a broad variety of small molecules (OH, 
NH, NO, SO2, ND3, CO, YbF, C7H5N, H2CO, LiH, CaH,...), but is currently limited 
to the production of molecules with a translational temperature around 10 mK. 

• Pairs of trapped ultracold atoms can be associated using laser flelds (photoassoci- 
ation, or PA) [3], or time- varying magnetic flelds (magneto- or Feshbach- associa- 
tion) [12]. Translational temperatures as low as a few tens of microKelvins can be 
reached, and even lower when quantum degeneracy is achieved. In contrast to the 
previous case, these results are up to now limited to alkali diatomic molecules which 
are most often formed in highly-excited vibrational levels, i.e. with high internal 
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energy. 

In the latter breakthrough occurred in 2008, when the possibihty to create 

ultracold bialkah molecules in the lowest vibrational level of their ground state or their 
lowest metastable triplet state, has been demonstrated. Caesium dimers have been 
created in their ground state f = level after a PA step followed by spontaneous 
emission and further vibrational pumping, using a sequence of shaped laser pulses |43j . 
Dipolar molecules, namely LiCs, have been observed in their absolute ground state 
rovibrational level f = 0, J = after PA and spontaneous emission [H]. Samples close 
to the degeneracy regime of ultracold molecules in the f = level of their ground state 
have been observed for Cs2 [IHl HH] , and KRb [57] , and in the f = level of their lowest 
triplet state for Rb2 [H], using the STIRAP technique (Stimulated Rapid Adiabatic 
Passage) to transfer the population from initial high-lying vibrational levels. 

2. Motivation of the present work 

The growing availability of ultracold samples of alkali diatomics brings the possibility 
to observe and study their interactions with neighboring ultracold atoms and molecules. 
In two independent - and almost simultaneous- experiments, inelastic rate constants 
for atom-molecule and molecule-molecule collisions have been extracted, using optically 
trapped Cs2 molecules created by PA and spontaneous emission |50]. The rate 
constants have been found independent from the initially populated rovibrational level. 
A resonant feature observed in the loss rate after collisions between trapped ultracold 
Cs2 molecules created by Feshbach association has been interpreted as Cs2-Cs2 bound 
states [51]. Inelastic atom-molecule collisions in a trapped sample of RbCs molecules 
have been studied for both RbCs-Cs and RbCs-Rb cases [52] . No systematic dependence 
on the internal state of the molecule has been probed within the experimental precision. 

Theoretical knowledge of the electronic structure of alkali triatomic systems 
is strongly needed for theoretical dynamical studies, which would support these 
experiments. There has been a significant research activity concerning the theoretical 
study of the structure and dynamical properties of homonuclear alkali-metal trimers. 
Non-additive effects in spin-polarized alkali-metal trimers were studied by Soldan et.al. 
[53] and the three-dimensional potential energy surfaces for the lowest quartet states 
were constructed for Lis [Ml ES], Naa [SH [57|, K3 [58] and Rbs ^9]. These were then 
used to study the corresponding spin-polarized reactive atom-dimer collisions at very 
low temperatures [551 EHl EOl EH [621 [63l EH EHl EHl [67] . General trends for the lowest 
quartet state of homonuclear alkali trimers - correlated to three spin-polarized alkali 
atoms - have been described by Hutson and Soldan [68] . One of the main characteristic 
of these systems is the importance of non-additive three-body forces, which are predicted 
to be large especially close to the equilibrium geometries [53|. Much less has been done 
in the case of heteronuclear alkali-metal trimers. One of us has actually started with 
a systematic study the lowest quartet states of Li2A mixed systems (with A=Na, K, 
Rb, Cs) [69]. He concluded that the single-reference coupled-cluster approach, which 
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has been successfully used for homonuclear alkali- metal trimers [SHI ES] , would be for 
various reasons very difficult to employ for calculations of the potential energy surfaces 
of heteronuclear alkali-metal trimers, and that an alternative approach to the problem 
should be sought. 

In the present paper, we present such an alternative approach. We extend our 
previous works on effective two-electron diatomic molecules like alkali dimers [TOl ED IZS] 
and alkali hydrides [73] to alkali trimers. We propose a preliminary study of the potential 
energy surfaces of the heaviest alkali trimer, Css, modeled as an effective three-electron 
molecule. Our approach is based on large-core effective core potentials (ECP) with core 
polarization potentials (GPP). We focused our work on two aspects: (i) the estimation of 
repulsion effects between the three Cs"*" ionic cores at short distances within the present 
ECP+CPP approach; (ii) the comparison of the results with the results obtained by the 
MOLPRO package for ab-initio calculations. 

3. Method of calculation 

3.1. Ah initio methods and basis sets 

The potential energies were calculated making use of the CIPSI package (Configuration 
Interaction by Perturbation of a multiconfigurational wave function Selected 
Iteratively) [74j. As in our previous studies on alkali dimers, the alkali atom A is 
described by an ^-dependent Effective Core Potential (ECP) for the ionic core A"*" [75|[76] 
including effective core polarization potential (CPP) terms [771 [78], and by a large set 
of uncontracted Gaussian functions for the valence electron. The atom is modeled as an 
effective one-electron system, permitting us to perform Full Configuration Interaction 
(FCI) calculations for any alkali dimer or trimer, considered as an effective two-electron 
or three-electron system, respectively. 

As the distance between different cores becomes smaller than the equilibrium 
distance of the system, it is well known that their electrostatic repulsion is not accounted 
for by this effective large-core potential. At the same time, the addition of the CPP 
term leads to an non-physical attractive behavior at short distances. For diatomic 
molecules, an empirical repulsion term is usually added to the calculation (see e.g. 
Refs. [79 l [801 ISTl [82]). If this short-range repulsive interaction in the triatomic molecule 
can be treated as the sum of pairwise interaction terms, computation of the potential 
surfaces would be greatly simplified. In order to check if such an approximation could 
be justified, we calculated the core-core repulsion at the Hartree-Fock level making 
use of the MOLPRO 2006.1 Quantum Chemistry Package [83]. For lithium and 
sodium, the all-electron correlation-consistent polarized core-valence quadruple-^ cc- 
pCVQZ basis sets [H] were used. For potassium, rubidium, and caesium, the small-core 
scalar relativistic effective core potentials ECPIOMDF, ECP28MDF, and ECP46MDF, 
respectively, together with the corresponding uncontracted valence basis sets [H3| were 
employed. 
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For the caesium dimer and trimer, the counterpoise-corrected dimer interaction 
energies were optimized using an algorithm implemented in MOLPRO. In this algorithm 
the interaction energies were calculated at the coupled-clusters level making use of a 
single-reference restricted open-shell variant [86] of the coupled cluster method [87] 
with single, double and non-iterative triple excitations [RCCSD(T)]. The small-core 
scalar relativistic effective core potentials ECP46MDF was employed together with its 
uncontracted valence basis set, which was augmented by one set of diffuse function in 
an even-tempered manner (aug-ECP46MDF). All electrons from the "outer-core" 5s5p 
orbitals were included in the RCSSD(T) calculations. The full counterpoise correction of 
Boys and Bernardi [88] was applied to all the interaction energies in order to compensate 
for the basis set superposition errors. 



3.2. Core Polarization Potential for a triatomic molecule 
The core polarization potential (CPP) is given by [77]: 

Vcpp = ^ fc ■ fc (1) 

c 

The electric field fc is the one produced at point Tc by all the other electrons and 
the ionic cores with static dipole polarizability ac- 



rc' 



where is the position vector from core c to electron Rcc' the one between cores c 
and d . The cutoff function /i(r, p) ensures the convergence of the integral over electronic 
coordinates. The Vcpp term contains a purely geometrical part Vnn depending only on 
the relative positions of the nuclei: 

f \ — ^ \ — ^ Rcc' ' Rcc" ry ry / r,\ 

For a diatomic molecule, Vnn is given by the familiar formula: Vnn = — -^-^^r^- 
For a triatomic system, we have: 



Zlai + Zla2 Zjai + Zfa^ Zfaa + Zja^ 

^-"-12 ^-"-13 ^-"-23 ^^"^ 

_ Z^Z-^ax R12 ■ R13 Z\Z-^a2 R12 ■ R23 _ Z\Z20iz R13 ■ R23 

-^12-^13 -^12-^23 -^13 -^23 

In those formulas, Zi and are, respectively, the net charge and static polarizability 
of ion core i. This term is independent of the ah initio calculation itself and is added a 
posteriori. 
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4. Volume effect for triatomic systems 

The part of the interaction between two nuclei other than the point charge repulsion has 
been dubbed "volume effect" by Jeung in his 1997 paper on alkah diatomics ^89j. Here, 
we are interested in this volume effect for alkali triatomic systems. More precisely, we 
investigate whether it can be reasonably approximated as a sum of two-body terms. 

The two-body repulsion energies Vcc are calculated for all the alkali homonuclear 
diatomic molecules. The three-body repulsion energies Vccc are calculated for 
particular geometries, namely for J^ooh and D^h symmetries, to provide a representative 
investigation. Calculations ar done at the Hartree-Fock level of theory starting with the 
orbitals of the free atoms. Those orbitals are kept frozen as the internuclear distance is 
shortened. In this respect, our calculations for the two-body term are similar to those 
of the "A" column from the tables of ref . [89] . This comparison is shown in Figure [U in 
logarithmic scale. Apart from a slight discrepancy in the case of Na2, there is a good 
agreement between both calculations, and it is seen that the two-body core repulsion 
term could be well fitted by an exponential form, as already stated in ref. [80] . 




Internuclear distance (Bohrs) 

Figure 1. Two-body core repulsion energies for allcali dimers. Blue : our work. Red 
dashed : Ref. [89]. Note the logarithmic scale used for the potential energy axis. 

Three-body core repulsion energies have been calculated for all homonuclear alkali 
trimers at linear symmetric and equilateral geometries according to the schemes of 
Figure [2l In Figure [3l we compare the three-body core repulsion term for linear 
geometries Vflf^{R) (the internuclear distance R being defined in Figure [2]), with the sum 
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of two-body term 2Vcc{R) + Kc(2-R). In the same manner, we compare in Figure H] the 
three-body repulsion term for equilateral geometries V^^{R) to 3V^'^\R). In both cases, 
it is seen that the three-body core repulsion term is well described by a pairwise additive 
lemma, up to the region where it becomes negligible. Therefore, just like for alkali 
dimers, the Vccc term can be added a posteriori to the potential surface calculations, 
whatever the chosen grid is for them. 




Figure 2. Linear (Doo/i) and equilateral (Dsh) geometries at which the three-body 
core repulsion term have been calculated. The Jacobi coordinates {r',R',d} used in 
the calculation of the potential energy surfaces of Css are also presented. 
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Internuclear distance (Bohrs) 



Figure 3. Additivity of the three-body core repulsion term for linear geometries. 
Blue : Three-body term V^^{R). Red dashed : sum of the relevant two-body terms 
2Vac{R) + Vcc{2R) (see text). 

5. Preliminary application to caesium trimer 

The inclusion of the three-body core repulsion term is considered in preliminary 
calculations concerning the state of Css which correlates to three spin-polarized 
caesium atoms. The used basis set is kept relatively small to maintain reasonable size for 
the FCl, reaching 152 532 Slater determinants in symmetry in which the calculation 
is performed. The basis set {dsSpid/ 4:s3p3d) is detailed in Table [H We use ^-dependent 
CPP's [78] which allow us to adjust the calculated atomic energies to the experimental 
ones for the 6s "^S, 6p^P and 5d atomic states of caesium. This is achieved by tuning 
the cutoff radii of Eq. [2] to the values displayed in Table [2j 

5.1. The triplet state of Cs2 

With these basis set and cut-off radii, we first calculate the lowest potential energy 
curve for Cs2 dissociating into Cs(6s)-|-Cs(6s) in the framework of our ECP-l-CC-FCl 
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^ 6 8 10 12 14 16 



Internuclear distance (Bohrs) 

Figure 4. Additivity of the three-body core repulsion term for equilateral geometries. 
Blue : Three-body term V^^^{R). Red dashed : sum of the relevant two-body terms 
3F(2)(^) (see text). 

method. We found an equilibrium internuclear distance Re = 11.82 oq and a well depth 
of De = 380.6 cm~^. These numbers can be compared to our optimized results obtained 
at the RCCSD(T)/aug-ECP46MDF level = 12.19 oq and De = 256.4 cm"! and the 
results from Ref. [53l [68] . who obtained Re = 12.44 ao and De = 246.8 cm~^ at the 
RCCSD(T)/ECP46MWB level of theory. The MOLPRO calculation yields a potential 
well which is deeper by is a more than 140 cm~^ than the one of the ECP-I-CC-FCI 
calculation. Apparently, the small-core ECP in combination with a rich valence basis 
set in the RCCSD(T) calculations provides more realistic description of the interaction 
energy than the large-core ECP in combination with a smaller valence basis set used in 
the FCI calculations (see also the discussion in Section [6]). The potential energy well 
of the triplet state results from the competition between the exchange energy at short 
distances and the attractive dispersion forces at long distances and is very sensitive to 
the quality of the basis. For this curve, the inclusion of the two-body core repulsion 
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Table 1. Gaussian basis set used on each caesium atom in the ECP+CPP-FCI 
approach. 
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2.6248 


p 


1.87 


d 


2.8111 



Table 2. ^-dependent cutoff radii pi (in units of oo) used in the core polarization 
potential. The polarization for Cs+ ionic core of 16.33 is taken from Ref [9D] . 

term does not change the characteristics of the well but ensures a correct exponential 
behavior at short internuclear distances (Figure [5]). 

5.2. Quartet state of Cs^ 

We then calculate the potential energy surface of the lowest quartet state of Css. 
Calculations are carried out in Jacobi coordinates {r', i?', 6} (depicted in Figure [2]) with 
the angle 6 fixed at 7r/2 in order to explore nuclear geometries, as a representative 
geometry. We have computed 984 points on a grid where 4.4 < r' < 40 Oq and 
4ao < -R' < 40 oq. The minimum of the surface is found at a D3/1 geometry with a 
bond length of R'^ = 10.72 oq and a well depth of De = 1518 cm~^ with the inclusion 
of the pairwise three-body core repulsion term, which ensures a realistic short-range 
repulsive wall. These numbers can be compared to our optimization results obtained 
at the RCCSD(T)/aug-ECP46MDF level R'^ = 11.18 ao and = 1208 cm'^ and to 
the results from Ref. [531 EH], who obtained R'^ = 11.33 ao and Dg = 1139 cm~^ at the 
RCCSD(T)/ECP46MWB level of theory. It is worth noting that as the minimum of 
the quartet potential well is located at shorter distances than the minimum of the Cs2 
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Figure 5. Calculated potential energy curves for the lowest state of Cs2. Black 
plain curve : ECP+CPP-FCI calculation with the basis detailed in Table [TJ Red 
dashed curve : obtained at the RCCSD(T)/aug-ECP46MDF level. Green chain dotted 
curve : ECP+CPP fuU CI result obtained with basis set "A" from Ref. [70] (see the 
discussion in Section[6]). Blue dotted curve : Multiparameter Morse Long Range fitting 
to spectroscopic data measured by Xie et al [91] . 



triplet state, the Vccc term indeed contributes to the depth of the well, and not only in 
the region of the repulsive wall. However, the magnitude of Vccc is still small, as the 
equilibrium internuclear distance is R'c = 10.7 cq and the well depth is De = 1522 cm~^ if 
we neglect it. As expected from the diatomic calculations above, we obtain a difference of 
about 300 cm^^ on the well depth of the quartet state of Css compared to the MOLPRO 
calculations. This confirms the limited quality of our basis set at the current level of 
our computations. 

5.3. Doublet states of Cs^ 

The three lowest potential surfaces for Csa are presented in Figure [61 The ^-82 and 
"^Ai states are the two components of a "^E' state at D3/1 symmetry subjected to 
Jahn- Teller effect. The left (resp. right) column shows the surfaces without (resp. 
with) the addition of the core repulsion term. Note the strongly unphysical attractive 
behavior in the region r' < 5ao, which is removed when Vccc is added. The global 
ground state ^i?2 is characterized by a well depth of 5437.1 cm~^ at the geometry 
{r' = 10.58 ao, R' = 7.36 Oq} with the inclusion of the core repulsion term. This changes 
to a well depth of 5489 cm~^ at {r' = 10.57 ao,R' = 7.3 ao} without the core repulsion 
term. Once again, the effect of Vccc is larger than in the dimer, as the equlibrium 
distance is shorter. The other component "^Ai is characterized by a well depth of 
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5258.2 cm ^ at {r' = 8.75ao)-R' = 9.28 ao} with the core repulsion term (5312 cm ^ 
at {r' = 8.66ao,-R' = 9.26 ao} without it). 




Radius r' of Jacobi (ao) 

Figure 6. Lowest three calculated surfaces for Cs3 in Jacobi coordinates with fixed 
6 = 7r/2. Left (resp. right) column : without (resp. with) the inclusion of the 
cores repulsion term Vccc- The red dashed curve locates the equilateral D3/1 geometry 
R' = ^r'. Contours are distant of 500 cm and the dashed contour is the energy of 
the dissociation Gs{^S) + Gs{^S) + Gs{^S). 

We present in Figures [7] and [8] the cut through the Ds/i geometry of the previous 
surfaces. We recall that in the D^h symmetry point group, the two doublet states 
and "^Ai discussed before correlates to the two components of a doubly degenerate '^E' 
state. We clearly see in these figures the abrupt attractive behavior which occurs at 
short distances and that the addition of the core repulsion term gives us back a realistic 
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Figure 7. Cut through the D3/1 geometry for the lowest quartet states of CS3. Plain 
black (resp. red dashed) curve : with (resp. without) the core repulsion term V^^^{R). 
Green chain dotted curve : quartet state on the RCCSD(T)/aug-ECP46MDF level of 
theory (see the discussion in section [51 




Figure 8. Cut through the D3/1 geometry for the lowest doublet states of CS3. Plain 
black (resp. red dashed) curve : with (resp. without) the core repulsion term V^^^{R). 
Both components of the ^E' state are represented. 
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6. Future prospects 

Preliminary potential energy surfaces for doublet and quartet states of the Csa molecule 
have been computed in the framework of a ECP+CPP-FCI quantum chemistry 
approach. The core repulsion interaction has been taken into account as an additional 
empirical correction as it is routinely done for diatomic molecules calculations. We 
showed that this term is well described by a sum of two-body terms. As mentioned 
before in ref. [HQ], it can be fitted by an exponential formula. Therefore, this term needs 
not to be calculated at each point of the three dimensional grid spanned by the degrees 
of freedom of a triatomic molecule on which the ah initio calculations are performed. It 
can easily be added a posteriori to cancel any non-physical behavior which might occur 
at short internuclear distances. 

To estimate the quality of this basis set, we have compared our calculations to high 
level ones which uses a different method than ours, with the MOLPRO package. Such 
a comparison between different methods is crucial to assess their accuracy and their 
consistency, as no spectroscopic data are available for the Css system, as well as for 
most of the triatomic alkali systems. The calculations presented here are exploratory 
in nature: the small basis set used (Table [1]) makes the calculation time short enough 
to generate over a thousand ah initio points. Figure [5] shows the limitation of this 
basis set where the well depth of the triplet state of Cs2 is overestimated by more 
than 100 cm~^ For the quartet state of Csa, our calculated well depth is consequently 
300 cm~^ lower than the RCCSD(T) calculations. If we employ the basis set previously 
used for the caesium dimer [TO] in the ECP-I-CPP-FCI method, this yields a very good 
agreement with the RCCSD(T) calculations, as the difference is now found at about 
20cm~^ between them (Figure [5]). The agreement is even more spectacular with the 
potential curve extracted from the spectroscopy of the lowest triplet state of Cs2 [91] . 

We are currently calculating a comprehensive set of ah initio points with the 
ECP-I-CPP-FCI with this extended basis for the Css molecule which will be the 
subject of a future paper. Such FCI calculations now involve more than 500,000 Slater 
determinants. In Figure [H we display the first calculated points for the Dooh symmetry 
of the lowest Css quartet state, compared to the MOLPRO calculation as described in 
Section [5l The figure suggests that the results are in good agreement with each other 
as the position of the minimum looks similar. The point calculated around 12.3 with 
the ECP-I-CCP full CI method is deeper than the minimum of the MOLPRO curve by 
about 40cm~^, while the difference between corresponding triplet potential curves in 
Figure [5] is seen to be about 20 cm~^ This suggests that pairwise additivity of the forces 
is a reasonably good approximation in this case. 
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Figure 9. Quartet "'E^ state of Css in Dooh symmetry. Green chain dotted 
curve : calculation at the RCCSD(T)/aug-ECP46MDF level of theory. Plus signs 
: preliminary calculations with the ECP+CPP-FCI approach using basis set "A" from 
Ref. [70]. 
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